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Abstract 

Ab initio calculations show that the coupling between antiferrodistortive(AFD) dis- 
tortions and magnetization in perovskite Bi2FeMn06 is prohibited to make mag- 
netization rotate as on-site Coulomb interaction([/) is larger than 2.7 eV, where 
anomalies in antiferromagnetic(AFM) vectors and band gap varying with on-site 
Coulomb interaction can be observed. This coupling is attributed to the antisym- 
metric Dzyaloshinskii-Moriya interaction (DMI) driven by the Cg-Cg states AFM in- 
teraction and charge redistribution with respect to different AFD distortions. 

Key words: Dzyaloshinskii-Moriya interaction(DMI);On-site Coulomb Interaction 
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1 Introduction 

Multiferroics are materials in which ferroic properties, e.g., magnetism 
and polar order coexist. Magnetic and ferroelectric ordering couple microscop- 
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ically or macroscopically to form the magnetic ferroelectrics. The couphng of 
the two ordering leads to the so-called magnetoelectric effect in which the mag- 
netization can be tuned by the external electric field, and vice versa [T]l2]l3] . 
Magnetic ferroelectrics have potential applications in information storage, 
actuators, sensors, and functional devices. Perovskite BiFeOa exhibits both 
weak ferromagnetism and ferroelectric characteristics, and has been studied 
extensively in recent days. |^|5|6|7f8] . The G-type antiferromagnetic (AFM) 
order of Fe magnetic moments exhibits a canting caused by the antisymmet- 
ric Dzyaloshinskii-Moriya interaction(DMI) under rhombohedral R3c space 
group. However, a spiral spin structure of AFM Fe sublattice rotates through 
the crystal with a long-wavelength period of 620A, and decrease the weak 
ferromagnetism further. In our previous work, we suggested that the effect of 
decreasing magnetism caused by the spiral spin structure can be suppressed 
by doping magnetic transitional metal ions in perovskite B sites[9]. We suggest 
perovskite Bi2FeMn06 is a good candidate to fulfill this requirement in that 
Mn and Fe have different magnetic moments and hence lower DMI to produce 
the rotation of AFM vectors. Comparing with BiFeOa, the new system should 
have an observed magnetization and reversal of magnetization macroscopi- 
cally. Except the ferroelectric displacement of Bi site driven by Bi-6s stereo- 
chemically active lone pair induced by the mixing between the (ns)^ ground 
state and a low-lying {nsY{npY excited state, there exists another structural 
distortion, the alternating sense of rotation of the oxygen octahedra along [1 
1 1] direction, which is known as antiferrodistortive(AFD) distortion [TOlITT] . 
Also in our former paper, it is shown that AFD distortion couples with the 
weak ferromagnetism due to the DMI which is decreasing with increasing on- 
site Coulomb iteraction([/) [12]. We aim to extend our study to perovskite 
Bi2FeMn06 to investigate the coupling between AFD distortion and weak 
ferromagnetism under DMI using ab initio calculations with considering the 
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spin-orbital(SO) coupling effect and the noncollinear spin configuration, while 
Mn were doped periodically along [1 1 1] direction in Fe site of BiFeOs .Does 
AFD distortion couple with magnetism in Bi2FeMn06 when Fe and Mn are 
arranged in G-type AFM ordering? What is the origin of coupling between the 
AFD distortion and the magnetism under DMI in Bi2FeMn06? What is the 
role of DMI in the coupling between AFD distortion and magnetization? How 
does Coulomb iteraction(?7) take effect on the DMI? In this paper we have 
proposed the mechanism of DMI in Bi2FeMn06, using first-principles calcula- 
tions based on density functional theory(DFT).This work can shed light on 
discovering the magnetoelectric coupling process in perovskite multiferroics. 

The remainder of this paper is organized as follows: In section 2, we 
presented the computational details of our calculations. We reported the cal- 
culated results and discussions in section 3. In section 4, we drew conclusion 
based on our calculation. 

2 Computational details 

Our calculations were performed within the local spin density approx- 
imation(LSDA) to DFT using the ABINIT package [T3|T4] . The ion-electron 
interaction was modeled by the projector augmented wave (PAW) potentials 
P^PB] with a uniform energy cutoff of 500 eV. Bi 5d, 6s, and 6p electrons, Fe 
4s, 4p,and 3d electrons, and O 2s and 2p electrons were considered as valence 
states. Two partial waves per I quantum number were used. The cutoff radii 
for the partial waves for Bi, Fe, and O were 2.5, 2.3, 1.1 a.u., respectively. 
6x6x6 Monkhorst-Pack sampling of the Brillouin zone were used for all cal- 
culations. We calculated the net magnetization per unit cell and the electronic 
properties within the LSDA+U method where the strong Coulomb repulsion 
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between localized d states has been considered by adding a Hubbard-like term 
to the effective potential [T7|18f 19] . The effective Hubbard parameter, the dif- 
ference between the Hubbard parameter U and the exchange interaction J 
[U — J), was changing in the range between and 5 eV for the Fe and Mn d 
states. For eV of (f/ — J), J was varying as 0,0.5, 0.8, and 1 eV, respectively. 
J remained 1 eV for other effective Hubbard values. Taking into account the 
SO interaction, we introduced the noncollinear spin configuration to construct 
the G-type AFM magnetic order with the AFM axis being along the x axis in 
Cartesian coordinates in our ab initio calculation. 

3 Results and discussion 

The initial lattice parameters are introduced as same as BiFeOa in Ref. [4] . 
Then cell shape and the atomic positions are fully relaxed, and the relaxed 
parameters are given in table 1. One can see that the lattice constant is de- 
creased with doping Mn in BiFeOs. The rhombohedral angle is also decreased 
comparing with BiFeOs. This shows that the lattice cell is compressed with 
substitution of Mn in Fe site while lattice shape is also changed with decreasing 
rhombohedral angle. 

For the AFD distortion, a rotational vector R has been introduced to 
describe the direction of the rotation of the oxygen octahedra[I2]. Rom* corre- 
sponds to the state in which the rotational vectors of two neighboring oxygen 
octahedra are deviating away, while Rj„ is pointing inward. The rotational 
angle is 10°in the Cartesian coordinates [12]. U and J are the Coulomb in- 
teraction and superexchange interaction parameter which were implemented 
in LSDA+U calculation as in ref. [12]. U is the amount of energy required to 
transfer one electron from one site to its nearest neighbor, while J indicates 
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the strength of the magnetic superexchange interaction. 

In table 2, we present the net magnetization per unit cell with respect 
to Rj„, and Hout in Cartesian coordinates for different U and J. For the sake 
of clarity, we only set varying value of J for U=0 eV. Again, one can see 
that J value have no effect on the resulting magnetization when U remains 
constant. The initial magnetic moment of Fe are arranged in ferrimagnetic 
order along x direction which is known as G-type AFM order. The resultant 
magnetic moment in x direction is due to the ferrimagnetic arrangement of 
Fe and Mn magnetization. The neighboring moment will interact under the 
DMI and deviate away from the original x direction to produce a resultant 
magnetic moment in y direction. When the rotational vector R is reversed, 
the resultant magnetic moment will be reversed too until the U is smaller than 
2.7 eV. This indicates the AFD and magnetic ordering is coupled and DMI 
is still active in this case. As U is greater than 2.7 eV, the coupling between 
these two ordering is prohibited. That is to say, the DMI is precluded in this 
case. It is worth mentioning that the component of magnetic moment in y 
direction for U =2.7 eV becomes zero. This indicates the critical value of U 
for depressing DMI is 2.7 eV in which the antisymmetric Fe and Mn magneitc 
moments stop interacting with each other. This critical value is smaller than 
that for BiFeOa where the DMI is turned off as U is approaching 2.9 eV|12j. 
This indicates that Fe-Mn antisymmetric interaction where interacting spin 
moment are (5/2)5' — (4/2)5" is weaker than Fe-Fe antisymmetric interaction 
where interacting spin moment are (5/2)5 — (5/2)5. We believe the weaker 
DMI and alternating spin moments in the new system would destroy the spiral 
spin structure through the crystal in comparison with the case in BiFeOa. 
Hereby the magnetization and reversal of magnetization can be improved in 
the new system. In order to make the reversal of magnetization clear, the 
resultant AFM vectors are illustrated in Fig.l with respect to U. AFM vectors 
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have three components corresponding to x, and z direction, respectively. 
It is apparent that there is an anomaly as V is approaching 2.7 eV which 
corresponds to the critical value of eliminating DMI. The magnetization is 
inverted as V is under the critical value which indicates the DMI is still active 
in this \J range. Meanwhile the DMI is inversely proportional to the \J value. 
The deviation of original magnetic moment away from x axis only occurs in 
xoy plane. That is why the reversal of magnetic moment can not be observed 
in z direction. 

The relationship between DMI and Coulomb on-site interaction can be 
understood theoretically as follows. The Hamiltonian of the new perovskite 
system reads, 

HBi2FeMn06 — ~2 ^ 3 Pei^Mnj^ Fei'^ Mnj~^ ^ ^ Fei,Mnj Fei^^ Mnj ■ 
<Fei,Mnj> <Fei,Mnj> 

(1) 

The first term is attributed to the symmetric superexchange, and the second 
term is the antisymmetric DMI contribution which only occurs when the in- 
version symmetry of the cations is broken. The constant D is related to the 
AFD displacement and the rotational vector. It functions like a 'DMI con- 
stant' which indicates the strength of the DMI. It reads by the second order 
perturbation in the case of one electron per ion 

J^PeMn = i^i/U)[bnn'{Fe-Mn)Cn'n{Mn-Fe)-Cnn'{Fe-Mn)bn'n{Mn-Fe)], 

(2) 

From this relation, it can be seen that the DMI constant is inversely pro- 
portional to the on-site Coulomb interaction and has nothing to do with the 
symmetric superexchange parameter. 

The new system should possess a good insulating property in terms of the 
band gap calculated, and hence a reversal states of ferroelectricity. The band 
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gap for metal oxides can be opened with adopting an appropriate U value. 
The corresponding value for BiFeOa is about 4.3 eV. We take 5 eV U value 
for Fe and Mn in Bi2FeMn06 to analyze the total density of states(DOS) in 
that the band gap can be increased slightly with increasing U value. The band 
gap corresponding to two rotational vectors are demonstrated in Fig.2. Gen- 
erally, the band gap corresponding to the two rotational vectors are increasing 
linearly with respect to U value, while there exist two anomalies to the two 
AFD states as U approaching 2.7 eV where the band gap are disappeared. 
These anomalies in the curve are consistent with the calculating results for 
AFM vectors shown in table 2 and Fig. 1, respectively. The anomalies show 
the disappearance of DMI occurs at the critical U value and is caused by the 
hopping of electrons through Fermi energy. Moreover, the band gap to Rout 
is more narrower than to Ri„ which is probably caused by the different lig- 
and field related with structure. In comparison with the total DOS without 
applying the on-site Coulomb interaction, the total DOS with U value of 5 
eV is shown in Fig. 3. The finite DOS in the vicinity of Fermi level is pushed 
away to form the band gap as U is equal to 5 eV. A 1.26 eV and 0.75 eV 
band gap is obtained for Hout and Rj„ rotational states, respectively. This 
indicates that these two AFD states both possess a semiconducting property. 
A ferroelectric reversal would be expected to these two reversal AFD states. 
However, Bi2FeMn06 remains metal property without applying U. Again this 
shows that the on-site Coulomb interaction is significant in opening the band 
gap for perovskite metal oxides. 

In order to shed light on the role of DMI between neighboring Fe and Mn 
ions, we report the orbital resolved density of states(ODOS) of the magnetic 
ions as U is equal to 2.6 eV which is exactly under the critical value for pre- 
cluding the antisymmetric DMI. The ODOS of Fe and Mn corresponding to 
Hout and Hin states are shown in Fig. 4. It is apparent that there are large 
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amount of doubly degenerate Cg states for Fe and Mn in the occupied bands 
around the vicinity of the Fermi level. The charge hybridization between Fe- 
3d and Mn-3d electrons is realized mainly by the Cg-Cg states interaction. It 
suggests that the Cg-Cg AFM interaction play an important role in the anti- 
symmetric DMI between Fe and Mn ions. Cg state is composed by dz^-r^ and 
dx^_y2 orbitals. It is worth mentioning that the finite ODOS in the vicinity of 
Fermi level is mainly attributed to the (i^2_^2 orbital. We suggest that Cg-Cg 
AFM interaction is exactly mediated by neighboring d^^^r'^ orbitals of Fe and 
Mn ions. The charge transformation between Fe and Mn ions, and hence the 
inversion of magnetization can be seen from Fig. 5. The hopping of electrons 
between neighboring Cg-Cg pairs occurs due to the DMI as the inversion sym- 
metry of cations is broken. Moreover the hopping of electrons couples with the 
rotation of the neighboring oxygen octahedra to produce the resultant magne- 
tization. The charge distribution in neighboring Cg-Cg states is polarized to a 
certain rotational vector, and the charge distribution is inversely polarized as 
the rotational vector is changing direction, leading to the reversal of resultant 
magnetization eventually. The charge distribution is supposed to distribute 
uniformly without applying the DMI. However, the homogenous distribution 
is destroyed by the DMI between neighboring magnetic ions to form a polar- 
ized status. We suggest this constitutes the microscopic mechanism of DMI 
in perovskite multiferroics. On the other hand, the AFM interaction in triply 
degenerate ^2^-^29 composed by d^y, dyz, and d^^ orbitals is relatively weak 
comparing with Cg-Cg AFM interaction. Again this shows that the latter is 
deeply involved in the antisymmetric DMI. 
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4 Conclusion 



The critical value of U for depressing DMI in Bi2FeMn06 is found to be 
2.7 eV which is smaller than that in BiFeOs. This indicates that the DMI in 
the new system functions weaker in comparison with that in BiFeOs . ODOS 
of Fe and Mn corresponding to different AFD displacements show that the 
microscopic mechanism of the DMI is originating from the hopping of electrons 
between neighboring Cg-Cg states, and the Cg-Cg AFM interaction couples with 
the rotation of the neighboring oxygen octahedra. 
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Table 1 

Calculated lattice constant a, rhombohedral angle a, volume V, and WyckofF 
parameters for Bi2FeMn06. The Wyckoff positions la(x,x,x) for the cations and 
3b(x,y,z) for the anions. 

BiaFeMnOe BiFeOs 



a(A) 5.030 5.459 

a(°) 59.02 60.36 

V(A3) 105.57 115.98 

Bi X -0.012/0.489 0.000 

Fe/Ti X 0.226(Fe)/0.727(Mn) 0.231 

O X 0.541 0.542 

y 0.950 0.943 

z 0.406 0.398 
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Table 2 

Magnetization per unit cell with respect to different value of U and J. 

U{eV) 0.5 0.8 1 2 

J{eV) 0.5 0.8 1 1 

-f^in R'out -E^in R'out f^in I^out R'in -R'out -f^in I^out 

Mr,{fiB) -2.9723 1.0683 -2.9723 1.0683 -2.9723 1.0683 -2.9723 1.0683 -2.9775 1.0041 

MyinB) -0.0079 0.0887 -0.0079 0.0887 -0.0079 0.0887 -0.0079 0.0887 -0.0100 0.0367 

M,{hb) 0.0004 0.0071 0.0004 0.0071 0.0004 0.0071 0.0004 0.0071 0.0021 0.0025 

U{eV) 3 3.5 3.6 3.7 3.8 
J{eV) 11111 

-f^in R-out I^in -f^out I^in I^out R'in -f^out R'in I^out 

M^ifiB) -2.9813 1.0159 -2.9843 1.0126 -2.9846 1.0108 -2.9851 1.0087 -2.9856 1.0118 

MyinB) -0.0050 0.0208 -0.0013 0.0194 -0.0004 0.0205 0.0000 0.0211 0.0002 0.0238 

M,{hb) 0.0021 -0.0008 0.0013 -0.0004 0.0011 0.0001 0.0009 0.0009 0.0003 0.0017 

U{eV) 3.9 4 5 6 
J{eV) 1111 

-f^in R-out R'in -f^out I^in -E^out R'in -R-out R'in I^out 

Mr^ifiB) -2.9862 1.0115 -2.9866 1.0115 -2.9919 1.0093 -2.9921 1.0074 

MyinB) 0.0007 0.0231 0.0011 0.0233 0.0090 0.0220 0.123 0.0144 

M^{hb) 0.0003 0.0025 0.0000 0.0025 0.0000 0.0018 -0.0003 0.0013 
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Figure captions: 

Fig.l Relationship between AFM vectors and U . 
Fig. 2 Band function of U. 

Fig.3 Total DOS for Rj„ and Hout witli U=5 eV and without applying U. 

Fig. 4 ODOS of Fe and Mn for Ri„ and R^^* with U=2.6 eV and without 
applying U. 

Fig. 5 Schematic diagram for rotation of magnetization and DMI between 
the neighboring Fe and Mn ions in Bi2FeMn06. The arrow denote the 
magnetic moment. 
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Fig. 1. Relationship between AFM vectors and U 
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Fig. 2. Band gap as a function of U. 
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3. Total DOS for Rj„ and Rout with [7=5 eV and without applying U. 
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Fig. 4. ODOS of Fe and Mn for Rj„ and Rout with [7=2.6 eV and without applying 
U. 




Fig. 5. Schematic diagram for rotation of magnetization and DMI between the neigh- 
boring Fe and Mn ions in Bi2FeMn06. The arrow denote the magnetic moment. 
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